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DISCLAIMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 USERS Group, in con- 
nection therewith. 



1620 USERS GROUP PROGRAM REVIEW AND EVALUATION 



^J^rogram No. Date_ 

Program Name: 



1. Does the abstract adequately describe what the program is and 

what it does? Yes No_ 

Comment 

2. Does the program do what the abstract says ? Yes Nc_ 

Comment 



3. Is the Description clear, understandable, and adequate? Yes No 

Comment 

4. Are the Operating Instructions understandable and in sufficient 
detail? 

Comment 

Are the Sense Switch options adequately described (if applicable)? 
Are the mnemonic labels identified or sufficiently 
understandable ? 
Comment 



5. Does the source program compile satisfactorily (if applicable)? Yes No 

Comment 

'o. Does the object program run satisfactorily? Yes No 

Comment 



Yes No 



Yes No 



Yes No 



7. Number of test cases run 

Are any restrictions as to data, size, range, etc. covered 

adequately in description? Yes No 

Comment 



8. Does the Program meet the minimal standards of the 1620 Users 

Group? Yes No 

Comment 

9. Please list any suggestions to improve the usefulness of the 
program. These will be passed on to the author for his consideration. 
Comment 



Please return to: 



Your Name 



Mr. Robert J. Robinson (PREP) 
Marquette University 
Computing Center 
1515 W. Wisconsin Avenue 
Milwaukee 3, Wisconsin 



Company 
Address 



User Group 
Code 



THIS REVIEW FORM IS PART OF THE 1620 USER GROUP ORG AMI Z ATTON 1 3 PROGRAM 



REVIEW AND EVALUATION PROCEDURE. 
TO PARTICIPATE IN THIS EVALUATION. 



NONMEMBER3 ARE CORDIALLY INVITED 



CORRELATION PROGRAM 
FOR SYMMETRIC AND ASYMMETRIC MATRICES 



Author: Mr. William E. Wright 

Southern Illinois University 

Data Processing Computer Center 

Carbondale, Illinois 



Deck Key 
Program Deck 



SOUTHERN ILLINOIS UNIVERSITY 
DATA PROCESSING AND COMPUTING CENTER 

CORRELATION PROGRAM 
FOR SYMMETRIC AND ASYMMETRIC MATRICES 
FOR THE IBM 1620 
PROGRAM ABSTRACT 



I, PURPOSE: This program, which will handle up to 999 variables, will compute 
all possible simple pairwise correlation coefficients among a set of variables 
or between two sets of variables. The following quantities are computed and 



II. LIMITATIONS: Variables may be from one to ten digits in length. There may- 
be as many as 999 variables and 99999 observations. Sums are accumulated as 
fixed-point numbers. Means and standard deviation coefficients are truncated to 
five decimals. Correlation coefficients are rounded to six decimals. 

HI. MATHEMATICAL METHODS: Standard Statistical Formula: 



memory is utilized. The correlation coefficient output is about 300 per minute 
for large problems. All other output is at full punch speed. Running speed is 
inversely proportional to the number of ij pairs, to the maximum variable 
size, and to the number of observations. As an example, about 2. 2 seconds 
are required to process observation in an intercor relation of 23 variables, of 
maximum size A. 

V. REMARKS: In problems involving a large number of variables, the data cards 
must be read in more than once. 



VI. EQUIPMENT SPECIFICATIONS: The following are required: card input/ 
output, 40K, indirect addressing, and move flag. 



onally punched: N, SX., 2X. 2 , S(X.X.), X., o\, r... 
3 r i i i j i i ij 




Programmed by 
William E. Wright 
October 1962 



SOUTHERN ILLINOIS UNIVERSITY 
DATA PROCESSING AND COMPUTING CENTER 

CORRELATION PROGRAM 
FOR SYMMETRIC AND ASYMMETRIC MATRICES 
FOR THE IBM 1620 



I. PURPOSE 

This program, which will handle up to 999 variables, will compute all 
possible, simple, pair -wise correlation coefficients among a set of variables 
or between two sets of variables, without missing data. The following 
quantities are computed and punched: N, 2X^, SX. 2 , S(X^X^), 1C, a\, r„. 

II. LIMI A T A TIONS 

A. Size of variables: 

Variables may be from one to ten digits in length. 

B. Number of variables: 

The maximum number of variables is 999. 

C. Number of observations: 

There may be as many as 99999 observations. 

D. Accuracy of results: 

Sums are accumulated as fixed-point numbers. Means and 
standard deviations are truncated to five decimals. Cor- 
relation coefficients are rounded to six decimals. 

in. MATHEMATICAL METHODS 




N-S(X.X.) - 2X.-2X. 

1 J 1 J 



r. . = 



1J N 2 <r. • <r. 



IV. STORAGE REQUIREMENTS 

The program, constants, work areas, add table, multiply table, and 
input area are stored in locations 00000 - 06725. The output area is location,' 



39920 - 39999. The word sizes are stored in locations succeeding 16725. 
The rest of the core is available for the summations. The storage required 
for inter cor relations equals V • (3 • SM + 10) +1/2 V • (V-l)- (2 SM + 5). The 
storage required for battery by battery correlations equals (VI + V2) . (3 • SM + 
10) + VI • V2 • (2 • SM +5). (V equals the number of variables, and 
equal the number of variables in battery 1 and battery 2, respectively, and 
SM equals the maximum variable size) . 

SPEED 

The speed is inversely proportional to the number of ij pairs, to the 
maximum variable size, and to the number of observations. Approximately 
2. 2 seconds are required to process one observation in an inter cor relation 
of 23 variables of maximum size 4. Sums are punched at a speed of approx- 
imately 125 cards per minute. The number of sums per card depends on the 
maximum variable size (sm)* For an SM of 4, approximately 5 pairs of means 
and standard deviations are punched per second. Correlation coefficients are 
punched at a speed of about 5 per second. 

INPUT CARD FORMS ^ 

A. Parameter card(s) 

(1) Inter c or relation 

Column numbers 

1-5 Numeric name of study 

6-8 Number of variables 

9-10 Size of variable 1 

11-12 Size of variable 2 

13-14 Size of variable 3 

etc. 

(2) Battery by battery correlation 

Column numbers 

1-8 Same as inter cor relation 

9-11 Number of first variable in battery 

1 (i. e. , if battery 1 begins with the 
fifth variable, 005 is the correct 
indication) . 



12-14 Number of variables in battery 1 

15-17 Number of first variable in battery 2 

18-20 Number of variables in battery 2 

21-22 Size of variable 1 

23-24 Size of variable 2 

25-26 Size of variable 3 

etc. 



If necessary, continue size indications on additional cards, beginning in 
column 1. A negative sign must be punched in column 1 of all parameter cards 
except the last one. (If there is only one parameter card, consider it the 
last one in this sense.) If all variable sizes are the same, place Program 
Switch 2 off and place the size in columns 9-10 for intercorrelations or 21-22 
for battery by battery correlations. 

B. Data cards 

Column numbers 

1-5 Ignored by the program but may be 

used for identification. 

6-7 The sequence number of the card in the 

observation , beginning with 00. Blanks 
may be substituted for zeroes. 

8- These are the values of the variables 

being correlated. Negative signs in- 
dicating negative values should be 
punched over the units position. Each 
variable must take up exactly the 
number of digits specified for it on the 
parameter card. If the variables for 
a single observation occupy more than 
one data card, do not split up the last 
variable on the card. Instead, begin 
anew in column 8 of the next card. 

The last data card in a study must be succeeded by a card with a 0-2-8 
punch in column 1. 

VII. OUTPUT CARD FORMS 

The first card punched contains the name of the study. The next card 
contains the number of observations. Sums and sums of squares are punched 



next, preceded by a blank card. Columns 1-3 of each card indicate "i" for 
the first variable whose sum is punched on that card. Sums and sums of 
squares are interlaced. The first statistic on each card is always a sum, 
and the last statistic on each card is always a sum of squares. Sums of 
cross products a re punched next, preceded by a blank card. Columns 1-3 of 
each card indicate M i M for the first cross product punched on that card. 
Column 4 is blank; columns 5-7 indicate 11 j". Means and standard deviations 
are punched next, preceded by a blank card. The form is the same as for 
sums and sums of squares. Correlation coefficients , preceded and suc- 
ceeded by a blank card, compose the last type of output. The form is the 
same as for sums of cross products. 

Preceding each statistic except number of observations are either a 
blank and a blank if positive, or a blank and a negative sign if negative. No 
statistic will be split onto two cards. All excess spaces at the end of each 
card will be blank. N is 5 digits; 2X^, 1?^, and o\ are SM + 5 digits. 2X^ 2 
and 2(X."X.) are 2 • SM + 5 digits. Correlation coefficients are 7 digits. 
Seven nines (9999999) will be punched for undefined correlation coefficients. 

VIII. EQUIPMENT SPECIFICATIONS 

The following features are required: card input/output, indirect ad- 
dressing, 40K, and move flag. The program may easily be altered for 20K 
or 60K. 

IX. TESTING 

This program has been tested using intercorrelations, battery by 
battery correlations, uniform variable sizes, non-uniform variable sizes, 
sizes ranging from 01 to 10, multiple studies, multiple pass studies, and 
multiple card observations. 

X. OPERATING INSTRUCTIONS 

Place program in read hopper. Follow this by parameter card(s), data 
cards, and a card with a 0-2-8 punch in column 1, for each study. Place 
blank cards in punch hopper. Set switches: 



Parity switch 

Input/output switch 
Overflow switch 



•program 



stop 



stop 



-5- 



Program Switch I 



of£- 



— symmetric matrix problem 
--battery by battery problem 
— uniform variable sizes 



on- 



Program Switch 2 



off- 



on- 



- -non-uniform variable sizes 



Program Switch 3 



off- 



— no punching of sums and sums 

of squares, 
--punching of sums and sums of 



on- 



Program Switch 4 



off- 



squares 

--no punching of sums of cross 



on- 



products 
--punching of sums of cross 
products 



Press Reset Key and Load Key. When program deck has loaded, press Start 
Key and Punch Start Key. The Reader Start Key must be pressed to read in 
the last cards when the read hopper has emptied. If "CARD ORDER" is typed 
out, the data cards are out of order. Beginning with the next-to-last card in 
the read stacker, order the data cards, place them in the read hopper, and 
press Start Key. If "REREAD DATA" is typed out for a study, replace the 
data cards for that study in the read hopper (without parameter cards), press 
Start Key and Reader Start Key, and continue as before. The message "END 
OF STUDY" will be typed out after the completed punching for each study. If 
for some reason, such as detection of an error, it is desired, to terminate 
processing, this may be done at any time by pressing the Instant Stop Key. To 
continue as for a new study, press Reset Key and Insert Key, type "4900402" 
and press Release Key and Start Key. 
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William E. Wright 
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nr 
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T M 1 


DS 
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n s 

lv o 
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DS 
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DS 
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1»80Z 
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n s 
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M T M F s 

in i inc. 


nr 
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A r! 


n s 
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XM1 


DS 


29 »TM1+4Z 




xcw 


DS 


27»TM2+1Z 




SP 


DS 


6»TM1-7Z 




ST 


DS 


6»TM1-13Z 




V 


DS 


3?A4+21Z 




VI 


DS 


3 ».A3+9Z 





00787 


00003 


V2 


DS 


3 9 A3+2 12 


00 7 99 


00003 


V3 


DS 


3 9 A 3 4 3 3 Z 


oozi 3 a 


00004 

V v v U 1- 


V4 


DS 


4»SSS152+34Z 


00927 


00003 


V5 


DS 


3 9 A5+9Z 


05997 


00002 


V6 


DS 


29 A30+75Z 


5 999 


00002 


V7 


DS 


2 » A30 + 77Z 


06 01 


00002 


V8 


DS 


2 » A30+79Z 


01753 


00002 


V9 


DS 


2 »CLEAR-1Z 


4. ft ft 


00004 


V 10 


DS 


4»B6+10Z 


049 04 


00004 


VI 1 


DS 


4 > B6 + 34Z 


L- U 9 5 7 


n n o o "3 

u u u u .5 


FBI 


DS 


3 » A5 + 2 1Z 


014 7 9 


n n n n q 

U U U U d 


p p 9 


1 J o 


3 , a 74.337 


n o r\ ^) !~\ 


U U U U 9 


W 


DS 


5>39920Z 


018 97 


2 


C M 


n ^ 


9 « A Q A '„ 9^7 


D 9 n ft Q 

U tut) 7 


000 2 


D 


DS 


2 ■» A10 + 35Z 


O O /i 7 "2, 

U U H f _0 


0^ 

u u u v..; J 


B I 


DS 


5 >SSS152 + 71Z 


U 9 O 9 


% 

<J U V v -J 


S 


DS 


5 > A 1 + 1 1 Z 


n r> ^ 7 7 


V., 1 V.. v . ■ .J 


si 


DS 


5 * A1+23Z 


n n ^ ft q 

l -' J 1 o > 


0000 5 


S2 


DS 


5 9 A 1 + 3 5 Z 


00A0 1 


C 5 


S3 


DS 


5 9 A 1 + 47Z 


7 9 


000 5 


S4 


DS 


5 » A4-13Z 


008 37 


00005 


S5 


DS 


5 » A3 + 7-1Z 


00 9 8 9 


00005 


S6 


DS 


5» A5+71Z 


^ 1 1 


0000 5 


S7 


DS 


5? A5+83Z 


013 13 


00005 


Q 


DS 


5 >ST0RE+47Z 


01487 


00003 


DDI 


DS 


3> A7+45Z 




00003 


DD2 


DS 


3»SUMS+9Z " 


01453 


00005 


DD3 


DS 


5 » A7 + 1 1Z 


042 89 


00005 


DD4 


DS 


5 » A25-33Z 


047 8 9 


00005 


BI I 


DS 


5 »B4+11Z 


02 3 32 


000C 5 


B IJ 


DS 


5» A12+66Z 


7 3 37 

v.- C ~) -J f 


000( 5 


E I I 


DS 


5»A12+71Z 


02 3 69 


0000 5 


E IJ 


DS 


5 » A9-1Z 


02 593 


00005 


EI 


DS 


5»A13+71Z 


030 25 


0000 1 


SW 1 


DS 


1» A18-13Z 


03 24 


00001 


SW2 


DS 


1 » A18-14Z 


n o q n i 


00001 


SW3 


DS 


1 > A12+.35Z 


03 3 45 


00005 


N 


DS 


5 »XI XI 2-37Z 


00449 


00005 


P 


DS 


5 * SSS152 + 47Z 


A 9 3 ^ 

^ D C. -J J 


00005 


PI 


DS 


5 » SIGN + 83Z 


06041 


00005 


P2 


DS 


5 9A31-69Z 


p, o 7 a o 


00005 


P I 


DS 


5 9 A36+59Z 


R 1 1 1 


' 00005 


PIJ 


DS 


5 ? A37 + 59Z 


03151 


00003 


I 


DS 


39 A.14 + 9Z 


06036 


00005 


T 


DS 


59A31-74Z 


06041 


00005 


K 


DS 


59A31-69Z 


C6235 


00005 


Kl 


DS 


5 »SI GN+83Z 


04813 


00005 


RS 


DS 


59B3+11Z 


06.235 


00005 


HUN 


DS 


5 9SIGN + 83Z 


00402 






DEND 


SSS152Z 
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SAMPLE INPUT 
BATTERY BY BATTERY CORRELATION 
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SAMPLE OUTPUT 
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005 

001 
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00 3 
00 3 
004 
00<+ 
00 5 
005 
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006 
00 7 
007 
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00 1 
004 
007 
001 
004 



0000273 71 
000028989 
000023790 

000024123 
000024123 
000015857 
000023581 



0000 1 87574797 
0000196191867 
0000168.154594 
0000173566019 
00001 73566019 
00000889^755 1 
0000254562739 



000035111 
000022839 
000022444 
000024589 
000024589 
000026375 



0000262265451 
0000160186907 
0000185028866 
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001 
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00 1 
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00 5 
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00 5 
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c - 
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0000173566019 
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0000078403628 
0000185803049 
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0000192473276 
0000165376178 
0000123682148 
0000070957738 
0000046452413 
0000078403628 
0000197128727 
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456780000 
^82460000 
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527500000 
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338178039 
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491780000 
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0000043737765 
0000048623030 
00 0.0 36082876 
0000036988725 
0000096288401 



0000185029097 

0000216104236 

0000141463679 

00001 105 023-87. 

0000117730677 

0000102420001 

0000094181316 

0000169919508 



177251159 
331547203 
390397087 
390397087 
5 3 54437Q1 



579780000 
4488 80000 



237147392 * * 
410566045 * * 



317140000 277916034 * 



!0l 



0101350 
1 -0 7017 76 



0548488 
0572395 



■ 0264495 
0360841 



0680578 
0801767 



0893323 
071.6651 



A 1 



SAMPLE INPUT 
JNTERGORRELATION 



1000102304040404040404040404040404040404040404010 20103 

SSS01001 11 122223333444455 556666777788889999 00001 1 li2 2 223§3 34444$ 5 5-566667 77?88«8 
SSS010199994111111111111 

SSS0200432 15432 65437464323 134679 136979463 23 1 13 134649797776463 13 197^797974434343 
SSS02 01 99995 99 999 9994 321 

SSS03002 13 13 164879 96467994497 8746 623 134 5497964 13 134 67 9 79 7443 13 169 7976 643 13 12 644 
SSS030 1999925555555§2131 

SSS0400979764641 3 136 797643 10 102 10 100 101054649 79764333 1000131667964348 7797431349 
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SSS0 5 0000 104043 738 3 99 393 82 8 2 8 171205052 71 75 88 8382 7278 69393 9357240240101057505 377 

SSS05 19999400 0000 00 0010 

Z 



SAMPLE OUTPUT 



10001 
00005 



001 
003 
005 
007 
009 

on 

013 

015 

017 
019 

021 

023 



000017370 

000027371 

000028989 

000023790 

000024123 
000015857 
000023581 

000013286 

000018290 
-000049995 

000000165 

000001665 



0000 120427832 

0000187574797 

0000196191867 

0000168154594 

0000173566019 
000008890755 1 
0000254562739 

0000176525956 

0000253671688 
0000499900005 

0000000012947 

0000001318347 



000021325 

000035111 

000022839 

000022444 

000024589 
000026375 
000012391 

000014864 

000004825 
000000017 

000000015 



0000102583949 

0000262265451 

0000160186907 

0000185028866 

0000197128727 
0000233642051 
0000104949331 

0000191642122 

0000135580859 
0O0O000000065 
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001 
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001 
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022 
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003 
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003 
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004 


00 5 


0000194868451 


004 
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000014836345 1 


004 


013 


0000199598480 



0000063663226 
0000067958645 
0000095317918 
0000111726774 
■0000173682630 
0000005622705 
0000 154999789 
0000085649645 
0000126327367 
0000088387330 
0000000071476 

000016385 8531 
0000136762371 
0000194794007 
0000105573286 
0000001 168365 
0000147485378 
0000192473276 
0000115840708 



0000117660004 
0000059867763 
0000086863602 
0000091 193386 
0000000049945 

0000118403406 
0000084593760 
0000103658069 
00001 19282368 
0000000736230 

0000151950009 
0000164019673 
0000098492839 
0000063526729 
0000000106215 
0000151742711 
0000120238714 
0000103146860 



0000104366207 * * 

0000055511421 * * 

0000050831036 * * 

0000125702788 * * 

00000005 57205 • * * 

0000076659123 * * 

0000104 5640 06 * * 

0000068241752 * * 

0000042666803 * * 

0000000066930 * * 

0000136941172 * * 

0000076706755 * * 

0000070275182 ■* *' 

-0000273682629 * * 

0000011789865 * * 

0000138792616 * n g * 

0000216104236 ' * * 

0000122681999 * * 



005 006 0000156992873 

005 CIO 0000165376178 

005 014 0000059477496 

005 018 0000020209904 

005 022 0000000084354 

006 007 0000132640493 
006 Oil 0000043737765 
006 015 000003979470? 
006 019 -0000228367161 

006 023 0000009635244 

007 008 0000173947713 
007 012 0000117730677 
007 016 0000062927050 

007 020 0000000092952 

008 009 0000169007784 
008 013 0000081906032 
008 017 0000017601613 

008 021 0000001239744 

009 010 0000078403628 
009 014 0000023665624 
009 018 -0000036686935 

009 022 0000000094391 

010 Oil 0000096288401 
010 015 0000119656405 
010 019 -0000245865411 

010 023 0000006480624 

011 012 0000106711389 
Oil 016 0000099093339 

011 020 0000000045338 

012 013 000023139618 1 
012 017 0000193715052 

012 021 0000001319296 

013 014 0000148121651 
013 018 0000139629590 

013 022 0000000115263 

014 015 000007647^737 
014 019 -0000123897609 

014 023 0000007551108 

015 016 0000153578438 

015 020 0000000021879 

016 C17 0000188549274 

016 021 0000001138115 

017 018 0000153786346 

017 022 0000000074865 

018 019 -0000048245175 

018 023 0000004819509 

019 020 -0000000169983 

020 021 0000000000649 

021 022 0000000001177 

022 023 0000000011877 



0000 130186629 
0000095776164 
0000 107868153 
■0000289861011 
0000009363294 
0000 125370525 
0000 1 1 0502387 
O000O6 1 577829 
0000000075045 

0000 163826015 
0000101769386 
^00003888342 1 
0000001246421 
0000046452413 
0000035328719 

-0000025204978 
OOOOOOOi i 2704 
0000036988725 
00000010 19333 
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